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Abstract 

A floating hemisphere under forced harmonic oscillation at very high and very low frequencies 
is considered. The problem is reduced to an elliptic one, that is, the Laplace operator in 
the exterior domain with standard Dirichlet and Neumann boundary conditions, so the flow 
problem is simplified to standard ones, with well known analytic solutions in some cases. 
The general procedure is based in the use of spherical harmonics and its derivation is based 
on a physics insight. The results can be used to test the accuracy achieved by numerical 
codes as, for example, by finite elements or boundary elements. 
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1 Classification (AMS/MSC) 

76Bxx Incompressible inviscid fluids 
76B07 Free-surface potential flows 
33C55 Spherical harmonics 
33F05 Numerical approximation 

2 Introduction 

A semi-numerical computation of the added mass for a floating hemisphere is shown, under 
an harmonic forced oscillation on the free surface of an irrotational, incompressible fluid, and 
linearized boundary conditions. Two standard problems are considered: the heave and the surge 
modes, that is, vertical and horizontal oscillations, respectively. For simplicity, our attention 
is restricted to a fluid of infinite depth. The added mass coefficients at very-high and very- 
low frequencies found in this work are in agreement with those given in literature by another 
strategies. 

As it is noted by Hulme Q , the hydrodynamic formulation of a floating hemisphere is analogous 
to the two-dimensional circular cylinder. The added mass coefficients are computed as = 
^kk/{py)i and the damping ones, as D'^^ = D^k/ {pVuj)^ where V = (2/3)7ri?^ is the hemisphere 
volume, p is fluid density, and uj is the circular frequency of the oscillation. The asymptotic 
values of these coefficients, for very slow and very high frequencies, can be obtained by analytical 
calculus, for instance, by a variable separation or image methods. For the surge/sway mode at 
very slow frequency, the boundary condition (j)^z = 0, where (j) is the velocity potential, is 
equivalent to a symmetry operation respect the plane z = and, then, corresponds to the 
solution of a sphere oscillating in an infinity medium. The added mass for the last case is half 
of the displaced volume, e.g. see [§, then, the surge/sway added mass coefficient is A'-^^ = 1/2, 
respect to the true displaced mass (2/3)7ri?^/>, where the half factor is due to the analytic 
prolongation. On the other hand, the asymptotic values of the added mass in heave mode 
are not too easy to obtain, and they may be obtained by a semi-numerical computation with 
spherical harmonics. 

The reason for doing this work is twofold. First, the method of solution adopted here is based by a 
physics insight rather a mathematical one. Second, the computation is near-exact, so the results 
can be used to test the accuracy achieved by numerical codes, as finite-element and boundary- 
element ones, adapted to wave-drag and seakeeping flow problems, e.g. see 0, ^ |ll|, |l2|, IC, 



Added mass computation for the added mass coefficients by Storti and D'EHa 



3 



3 An oscillating hemisphere 

An oscillating hemisphere in a forced motion is considered. The flat face of the hemisphere is 
on the free surface of an irrotational and incompressible fluid, without a mean flow, where the 
fluid depth is assumed as infinity. The upward direction is z and the hydrostatic equilibrium 
plane is z = 0. Due to the symmetry, a spherical coordinate system is chosen as 

' z = r cos B ; 

< a: = r sin cos ; (1) 
,y = r sin 9 sin 4> 

At very high frequencies, the free surface boundary condition shrinks to the simple form (f) = 0. 
But, by anti-symmetry respects to the plane z = 0, the flow problem is equivalent to solve the 
modifled one 

' = inn' ■ 

< (f)^n = cos e at ; (2) 
> = at r[) ; 

where (f) is the velocity potential, A is the Laplace operator, Fg and Fq are the extended 
flow domain, extended hemisphere surface and extended free surface, respectively, through the 
reflection plane z = 0. As the free surface boundary condition for very high frequencies is 
(f> = 0, then, its right hand side term has been extended in an anti-symmetric way. That is, the 
extended problem at very high frequencies is the same that a sphere in infinite medium. On the 
other hand, the free surface boundary condition for very slow frequencies is (f)^n = 0; and for the 
extended problem, its right hand side must be extended in a symmetrical way and, then, 

' A(j) = in J)' ; 

< = I cos 6*1 at F^ ; (3) 
. ^ = at F(, ; 

where, due the module on | cos 6\, the slow frequency radiation problem does not have, in general, 
a closed solution and, then, it must be found with another resources like spherical harmonics, 
as it is considered in this work. 



4 The very-high and very-low frequencies limits 

The free-surface boundary condition in the limits of very-low and very frequencies 

(t>n = — </> ; (4) 
9 

is reduced to the homogeneous Neumann and Dirichlet boundary conditions, respectively, where 
g is the gravity acceleration. Also, the radiation boundary condition at infinity imposes that 
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Figure 1: Geometrical description of a seakeeping-like flow problem. 

the velocity potential (p tends to zero, so the corresponding Partial Differential Equation (PDE) 
system is 

f A(/) = in r2 ; 
6 „ = iioh at ; 

= at Fi? for low frequencies; (5) 
6 = at FiT' for high frequencies; 
^(f)\ for [x| ^ oo ; 
where the load h is the normal displacement of the mode under consideration. It can be seen 
that, under these conditions, the flow problem is transformed in a standard elliptic one, whose 
solution is real valued (in reality it is an imaginary one but it can be transformed by means of a 
simple re-definition). It is assumed that the load h is real (that is, the body motion is in phase). 
The added mass for the mode motion is found from 

ajj = ~~2 dTiucjy^n ^ ~ dTil^^n ] (6) 

where now ^p = —(i/uj)(f), so in the limits u; ^ and uj — > oo, the function ip is real. 
By symmetry, the Eqns (^ c-d) can be reproduced extending the flow problem to z > 0, by 
means of a mirror body image and extending the load h in an appropriate way. For instance, 
the homogeneous Neumann boundary condition is obtained extending the load in a symmetrical 
way with respect to z = 0, that is, 

Hx,y,z) = +h{z,y, -z) ; (7) 
while the Dirichlet boundary condition can be obtained extending in a skew-symmetrical way 

Hx,y,z) = +h(z,y, -z) ; (8) 
For example, the boundary boundary condition for the hemisphere in the heave-mode is 



tp^n = COS ; 



(9) 
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Figure 2: Extension of tlie load h for the heave- mode at very-low and very- high frequencies. 
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Figure 3: Extension of the load h = cos0 for the surge-mode at very-low frequencies. 



(10) 



where (r, 99, 0) is the spherical coordinate system with origin in the center of the hemisphere 
such as 

z = r cos i-p 
X = r sin (p cos 6 
y = r sin 93 sin 9 

here the Hildebrand's convention [Q] is used. Then, the limits at very-low and very-high fre- 
quencies for the hemisphere can be computed from the sphere ones, but with the loads 

{V',n = I cos 93 1 for low frequencies ; 
'(/',n = cos (f for high frequencies ; 

see figure |2|. Similarly, for the surge mode (oscillation along the x-axis) the equivalent load is 



11) 



h {ifO) = cos 6 ; 

and the extensions for very-low and very-high frequencies are 

{ip^n = cos 9 for low frequencies ; 

ip^n = sign (cos 93) cos 9 for high frequencies ; 



(12) 



(13) 
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z high frequency 
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h ((p,9) = sign [cos (pi cos 9 



Figure 4: Extension of the load h = sign (cos ip) cos 9 for surge-mode at very-high frequencies. 

5 Solution of the flow problems 



The previous flow problems can be solved, for instance, by an analytical way or by series. The 
heave solution at very-high frequencies and the surge one at very-low frequencies are the same 
of a sphere in an infinity medium and uniform velocity, so the additional mass is a half of the 
displaced mass, that is, 

J 033(0; ^ oo) = ; 

\ aii(a; ^ 0) = fpi?3 . 
In the other two cases it is necessary to make an expansion of the sources by means of spherical 
harmonics. 



(14) 



6 Spherical harmonics 



The solution of the exterior potential problem 

( A'(p = for r > 1 ; 
l^ = /(<^,0) atr = l; 

where ip = ip{<f, 9), it can be solved expanding the function f{ip, 9) in terms of the harmonics 



(15) 



00 n 



f{v>,G) = X] "^"O Pn{coscp) + X! X! [^""1 cos(m6') -I- bnm sin(m6')] P^{cosip) ; (16) 

n=0 n=0 m=l 



where 



2n + l 
Ait 



[ dTf{^,9)Pn{cosip); 

Jr=l 



(17) 



anm = ^^ i"^""!; / dr/(<^,^)P-(cos<^)cosme; 
ZTT (n + m)! Jr=i 



(18) 



6nm = ^^j^-^ / dr/(^,e)Pr(cos^)sinme; 
ZTT [n + m)l Jr=i 



(19) 
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where dT = sin (p dO dip is the soUd angle differential in spherical coordinates. Once this expan- 
sion is computed, the exterior potential can be written as 

oo 

V'(r,^,^) = ^y„(^,^)r-("+i); (20) 

n=0 

where 

n 

Yn{}P, 0) = anoPn{cos p) + ^ [anm cos{m9) + bnm sin(m6l)] -P™(cos if) ; (21) 

m=l 

now, the Neumann problem can be solved taking derivatives with respect to r and evaluating 
at r = 1, obtaining the expression 

oo 

h{ip,0) = -J2in + l)Yn{ip,0) ; (22) 

n=0 

from which analogous relations are obtained 

2n + 1 f 

^nO = -n — r-jT / dV h{ip,9) Pnicosp) ; (23) 

47r(n + 1) Jr=i 

anm = I K^^ C^(cos cos{m9) ; (24) 

2n{n + Ij (n + m)\ Jr=i 

_ 2n+l (n-m)! 



Once obtained the coefficients of the expansion, the additional mass is obtained from 

Ajj = f dr vV'.r ; (26) 

Jr=l 

= / difl al,[Pn{cosip)f + ^ (aL + 6L)[^r(cos^)]2 ; (27) 
where the orthogonality property of the spherical harmonics was taken into account, and 

Ajj = [ dr 4} V,r (28) 

Jr=l 

= /'^/^|«no[^n(cos^)]2+ f^(aL + &D[C^(cos^)]4 . ^29) 



m=l 



and using the properties of the Legendre polynomials 



oo 



, 2n + 1 

n=l 



T (n + m)\ 

m=l ^ ' 

7 Hemisphere in heave at very-low frequencies 



(30) 



In this case the load h{(p, 9) = \ cos (p\, so 

_ 2n + 1 
''""47r(n + l)y_i 



j_^dii\ii\Pn{lj); (31) 
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as tlie Pk are even (odd) for k even (odd), only remains the even terms and then 

^" = o \ ^ [ fJ- dfJ. PnifJ') for n even. (32) 
27r(n + 1) Jo 

For computing the integral, the P„ terms are generating in a recursive way from Pq = 1, Pi = fi 
and P2, Pn are obtained solving 

(n + 1) Pn+i - (2n + 1) ^ Pnifi) + n P„_i(/u) = . (33) 

The coefficients of the polynomials fiPn{fi) are obtained from the Pn ones, and the integral is 
made in a direct way. The final result is 

ass = 1-7403 p ; (34) 

corresponding to 033 = 0.83093, that is, the non-dimensional coefficient respects to the hemi- 
sphere mass 2/3TrpR^. 

8 Hemisphere in heave at very-high frequencies 

In this case, the factor cos 9 made that the only non-null coefficients are the a^i ones. For 
obtaining these, an integral from ^ = 1 to — 1 must be done, with a function that contents the 
Pn ones. These terms having a factor ^1 — /i^, so it is convenient to compute the integral in a 
numerically way. Our final numerical result is 

an (w ^ 00) = 0.29806... (0.14231...); (35) 

where the value between parentheses is referred, always, to the non-dimensional value respect 
to displaced mass. 

In brief, our estimates for the asymptotic values at very-low and very-high frequencies (subindex 
1 for the surge mode and 3 for the heave one), of the non-dimensional coefficients respects to 
the hemisphere mass 2/37rpi?3 (denoted with primes) are 

r A[i 0.5 for Ka^O ; 



A[i 0.14231 for Ka^ 00 
0.83093 for Ka^O ; 



(36) 



, ^433 0.5 for Ka 00 . 

On the other hand, some literature values found for the surge/sway mode, e.g. see Sierevogel 
Q, Prins 1^ (where only the intervals [0.25, 1.50] and [0.6, 1.5] are considered, respectively, so 
the extrapolations are rather doubtful), and for the heave one, e.g. see Korsmeyer and Liapis 



1 13(1, are 

r A[i 0.50 for ; 
A',, 0.25 for Ka^ 00 ; 

(37) 
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and the Hulme's Q ones are 

; 

oo ; 

oo ; 

Korsmeyer used a panel method with Fourier transform and complex impedance extended to 
very slow frequencies, while the Hulme's numerical results are obtained by spherical harmonics. 
The Sierevogel, Prins and Liapis results are obtained with a panel method and Kelvin source. 
In general, the concordance between our estimates and the literature ones is good, except in our 
surge mode coefficient A'n, which is about a half. 

9 Conclusions 

We have considered a semi-numerical computation by spherical harmonics for added mass co- 
efficients of an oscillating hemi-sphere at very low and very high frequencies. The method is 
based by a physics insight and the computation is near-exact, so the results can be used to test 
the accuracy achieved by numerical codes, as finite-element and boundary-element ones. 
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